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Abstract. In view of recent claims that asteroseismology could supply invaluable insights into the crystallization 
process occurring in the interiors of massive white dwarf stars, we present in this work new pulsational calcula- 
tions for improved carbon-oxygen DA white dwarf models suitable for the study of massive ZZ Ceti stars. The 
background models employed in this study, presented in detail in a recent paper by Althaus et al. (2003), are 
the result of the complete evolution of massive white dwarf progenitors from the zero-age main sequence through 
the Asymptotic Giant Branch (AGB) and mass loss phases to the white dwarf regime. Abundance changes are 
accounted for by means of a full coupling between nuclear evolution and time-dependent mixing due to convection, 
salt fingers, and diffusive overshoot. In addition, time-dependent element diffusion for multicomponent gases has 
been considered during the white dwarf evolution. Crystallization and chemical rehomogenization due to phase 
separation upon crystallization in the core of our models have been fully considered. The efi^ects of crystallization 
on the period spectrum of these massive white dwarf models are assessed by means of a detailed pulsational 
analysis of linear, nonradial, adiabatic gravity modes. To properly account for the effects of the presence of a solid 
phase in the models we employ special conditions on the oscillation eigenfunctions at the solid-liquid interface. We 
find that the theoretical pulsation spectrum is strongly modified when crystallization is considered, in particular 
concerning the mode trapping properties of the equilibrium models. We show that the strong mode trapping seen 
in the models with overshooting can be reproduced by means of a simple analytical model. We also discuss at 
some length the implications of our study for BPM 37093, the most massive ZZ Ceti star presently known. In 
particular, we attempt to place constraints on the physical processes occurring prior to the formation of this white 
dwarf. We find that if BPM 37093 has a stellar mass of ~ 1.00 Mq its observed spectrum could bear the signature 
of overshoot episodes during the helium core burning. 
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1. Introduction 

ZZ Ceti (or DAV) stars are cool, hydrogen-rich pulsat- 
ing white dwarf stars belonging currently to the most 
extensively studied family of degenerate pulsators — see 
Gautschy & Saio (1995, 1996) for a review. Their pulsat- 
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ing nature is evident from periodic brightness variations'^ 
caused by spheroidal, non-radial (7(gravity)-modes of low 
degree {£ < 2). Excitation of 100-1200 s oscillation modes 
is likely to be due to the action of the so called convective 
driving mechanism (Brickhill 1991, Goldreich & Wu 1999), 
although the early works of Dolez & Vauclair (1981) and 
Winget et al. (1982) postulated the k — j mechanism as 



^ Also fiuctuations in radial velocity have recently been mea- 
sured in some ZZ Ceti stars, see van Kerkwijk et al. (2000) and 
Kotak et al. (2002). 
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Fig. 1. Abundance by mass of ^H, ^He, ^^C and ^^O (up- 
per panel) and the run of the Brunt- Vaisala frequency 
(lower panels) as a function of the outer mass fraction q 
for the O.936-M0 white dwarf remnant of 7.5-Mq evolu- 
tion (sequence NOV) near the ZZ Ceti instability strip. 
The model is characterized by log(L/L0) = —2.92 and 
logrefi = 4.06. 

responsible for driving. ZZ Ceti stars are well known to 
pulsate in a "pure" instability strip at effective tempera- 
tures (Toff) between about 11000 and 12400 K. Relevant 
studies devoted to exploring the pulsational properties of 
ZZ Ceti stars are those of Brassard et al. (1991b; 1992ab), 
Gautschy et al. (1996), Bradley (1996, 1998a, 2001) and 
Corsico et al. (2001, 2002) amongst others. 

Over the last years, pulsation studies of ZZ Ceti stars 
— as well as of DBV, DOV and PNNV stars, the other 
categories of pulsating white dwarfs — through astero- 
seismology have become a valuable technique for sounding 
the white dwarf interiors and evolution. Indeed, asteroseis- 
mological inferences have provided independent valuable 
constraints to fundamental quantities such as core compo- 
sition, outer layer chemical stratification and stellar mass 
(Pfciffcr et al. 1996; Bradley 1998b, 2001). 

An outstanding application of white dwarf asteroseis- 
mology that has drawn the attention of researchers is re- 
lated to the possibility of placing observational constraints 
on the crystallization process in the very dense interiors of 
white dwarfs. Crystallization in the core of white dwarfs 
was predicted theoretically about 40 years ago (Van Horn 
1968; see also Kirzhnits 1960; Abrikosov 1960 and Salpeter 
1961), but it was not until recently that researchers have 
been able to peer into the crystallized structure of white 
dwarfs. This has been prompted by the discovery of pul- 
sations in the star BPM 37093 (Kanaan et al. 1992), a 
massive ZZ Ceti star which should be largely crystal- 
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Fig. 2. Same as Fig. [T]but for the O.94-M0 white dwarf 
remnant of 6-M0 evolution (sequence OV). Core over- 
shooting leads to a inner chemical profile quite different 
from the situation in which no extra mixing beyond the 
fully convective core is considered (Fig. In particu- 
lar note the almost-discontinuity in the chemical profile 
at log(l - Mr/M,) w -0.7. 

lized (Winget el al. 1997). In a very detailed investigation, 
Montgomery & Winget (1999) have explored the effects 
of crystallization on the period pattern of massive white 
dwarf stellar models. Amongst others findings, these au- 
thors concluded that some periods are notoriously sensi- 
tive to changes in the crystallized mass fraction of their 
models. They concluded that pulsating white dwarfs are, 
in principle, very promising objects to place constraints 
on the crystallization processes in stars. 

Very recently, Metcalfe, Montgomery, & Kanaan 
(2004) have performed asteroseismological fits to the pul- 
sation periods of the star BPM 37093. In their models, 
the crystallized mass fraction was treated as a free param- 
eter, with the goal of empirically determining the degree 
of crystallization for a given stellar mass and core compo- 
sition. Their preliminary results are that BPM 37093 has 
a crystallized mass fraction greater than 50%, with values 
possibly as high as 90%, although they indicated that an 
exploration of a finer grid of models (in both stellar mass 
and crystallized mass fraction) will be necessary in order 
to make conclusive statements. 

Crystallization has consequences for the car- 
bon/oxygen distribution within the core of a white 
dwarf. As a matter of fact, solid theoretical evidence 
suggests that, when crystallization occurs, the oxygen 
content in the solid phase is enhanced relative to that in 
the original fluid phase (Stevenson 1980; Garcia-Berro et 
al. 1988; Ichimaru et al. 1988; Segretain & Chabrier 1993). 




Fig. 3. Evolution of the periods for t — 2 modes as a Fig. 4. Same as Fig. 13 but for the OV model-sequence, 

function of the effective temperature, corresponding to the Thick lines emphasize some examples of modes experienc- 

NOV model-sequence. The effects of crystallization have ing bumping/avoided crossing phenomena, 
been neglected in computing the eigenmodes. 



As a result, the carbon content in the fluid surrounding 
the solid core is enhanced. Since carbon is lighter than 
oxygen, these fluid layers are Rayleigh- Taylor unstable, 
and the ensuing convective mixing will redistribute the 
abundances and lead to flat profiles in a region whose 
size depends on the initial composition profile and on 
the degree of chemical enhancement produced during 
the crystallization process (Isern et al. f997; Salaris 
et al. 1997; Montgomery et al. 1999). These mixing 
episodes could be particularly relevant in the context 
of carbon/oxygen cores with varying chemical profiles 
since they, in principle, could be smoothed out by such 
chemical rehomogenization. To the best of our knowledge, 
pulsational studies of massive white dwarfs taking into 
account the chemical rehomogenization induced by phase 
separation have not been performed; this effect was not 
considered in the pulsational calculations of Montgomery 
& Winget (1999), although it was fully taken into account 
in the calculations of phase separation and cooling delay 
in white dwarfs by Salaris et al. (1997), Montgomery et 
al. (1999), and Salaris et al. (2000). In view of these con- 
siderations, one of the main motivations for the present 
paper concerns the full exploration of the pulsational 
properties of massive white dwarfs that takes into account 
the effects of phase separation in carbon/oxygen cores. 

Recently, Althaus et al. (2003) — hereinafter Paper I 
— have presented new stellar models for massive ZZ Ceti 
stars, based on detailed evolutionary calculations that ac- 
count for the complete evolution of the white dwarf pro- 
genitor. In particular, such calculations include diffusive 



mechanical overshooting during the core burning phases 
of the progenitor star. In addition, the evolution of the 
chemical abundance distribution due to time-dependent 
element diffusion during the whole white dwarf regime has 
been considered in Paper I. In that work we concluded 
that, as a result of the smoothness of the chemical profile 
caused by diffusion processes, the mode trapping due to 
outer chemical interfaces is notably diminished, irrespec- 
tive of the occurrence of core overshooting. Instead, the 
theoretical pulsational spectrum is characterized by the 
presence of pronounced non-uniformity in the spacing of 
consecutive periods in models with core overshooting, at 
variance with the situation in which this mixing process 
is neglected. In particular, we found that the pulsational 
spectrum in models with overshooting are dominated by 
the presence of "core trapped" modes, characterized by 
relatively high values of the oscillation kinetic energy and 
strong minima in the period spacing diagrams. 

The main conclusion drawn in Paper I is that the pul- 
sational properties of massive ZZ Ceti stars become very 
sensitive to the occurrence of extra mixing episodes that 
take place beyond the formally convective core during the 
central helium burning, for instance core overshoot. Such 
mixing episodes give rise to a sharp variation of the core 
chemical composition that leave noticeable signatures on 
the theoretical period spectrum of pulsating white dwarfs. 
However, the fact that the eigenfunctions of g-modes are 
expected to have very low amplitudes in the solid core 
— because of the non-zero shear modulus of the solid 
(Montgomery & Winget 1999) — as well as the mixing 
episodes induced by crystallization were two effects not 
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Fig. 5. Same as Fig. O but for the case in which the 
presence of a solid core has been considered in the com- 
putation of the eigenmodes. Chemical rehomogenization 
upon crystallization has not been taken into account in 
the stellar model. The upper scale measures the crystal- 
lized mass fraction. 



addressed in Paper L Therefore, an initial question to be 
answered in this paper is whether the chemical disconti- 
nuity caused by core overshoot could be wiped out by the 
chemical rehomogenization, and more importantly, what 
implications this would have for the pulsational proper- 
ties of massive ZZ Ceti stars. In this connection, another 
aim of our work is to place constraints on the stellar mass 
and Toff values at which we should expect a pulsational 
pattern without any signature of core overshoot. 

We believe that a re-examination of the pulsational 
properties of massive white dwarf stars deserves to be 
done in the frame of our new massive ZZ Ceti model stars. 
Specifically, our aim is to explore the pulsational proper- 
ties of the white dwarf models presented in Paper I tak- 
ing crystallization self-consistently into account. Section 
§2 contains a brief description of the main physical inputs 
of the models. In §3 we explore the response of the pul- 
sational pattern to the presence of crystallization. Section 
§4 is devoted to discussing the implications of our results 
for observational expectations and, finally, in §5 we make 
some concluding remarks. 

2. Input physics and evolutionary sequences 

Our massive, carbon/oxygen-core white dwarf models ap- 
propriate for ZZ Ceti stars have been obtained with the 
LPCODE evolutionary code described at length in Paper I 
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Fig. 6. Same as Fig. |31 but for the OV model-sequence. 



and references therein. The code is based on an up-to-date 
and detailed physical description, particularly regarding 
the processes responsible for the chemical changes. In this 
section, we summarize the input physics of our stellar code 
and the main characteristics of our evolutionary models. 

LPCODE uses OPAL radiative opacities — includ- 
ing carbon- and oxygen-rich compositions — for arbitrary 
metallicity from Iglesias & Rogers (1996) and molecular 
opacities from Alexander & Ferguson (1994). The equa- 
tion of state includes partial ionization, radiation pres- 
sure, ionic contributions, partially degenerate electrons 
and Coulomb interactions. For the white dwarf regime, 
we include an updated version of the equation of state of 
Magni & Mazzitelli (1979). Neutrino emission rates and 
high-density conductive opacities are taken from the works 
of Itoh and collaborators — see Althaus et al. (2002). 
A total of 34 thermonuclear reaction rates and 16 iso- 
topes characterizes our nuclear network that describes the 
hydrogen — proton-proton chain and CNO bi-cycle — 
and helium burning, and carbon ignition. Nuclear reac- 
tion rates are taken from Caughlan & Fowler (1988) and 
Angulo et al. (1999). 

As for the abundance changes, we consider a time- 
dependent scheme for the simultaneous treatment of 
chemical changes caused by nuclear burning and mixing 
processes. Specifically, the changes in the abundances for 
all chemical elements are described by the set of equations 
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with Y being the vector containing the number fraction 
of all considered nuclear species. Here, mixing due to con- 
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Fig. 7. The forward period spacing {APk — Pk+i — Pk) in 
terms of the periods for i — 2 modes corresponding to the 
NOV model-sequence. Each panel shows the predictions 
for different percentages of the crystallized mass fraction. 
Chemical rehomogenization was not considered. 



vection, salt finger and overshoot is treated as a diffusion 
process which is described by the second term of Eq. iQJ . 
The efficiency of convective and salt-finger mixing is de- 
scribed by appropriate diffusion coefficients, D, which are 
specified by the treatment of convection. In particular, 
we considered the extended mixing length theory of con- 
vection for fluids with composition gradients developed 
by Grossman et al. (1993) in its local approximation as 
given by Grossman & Taam (1996). This treatment applies 
in convective, semiconvective^ and salt finger instability 
regimes. The first term of Eq. gives the abundance 
changes due to thermonuclear reactions, changes which 
are fully coupled to mixing processes. This term is lin- 
earized following the implicit scheme of Arnett & Truran 
(1969). For details concerning the numerical scheme em- 
ployed to solve Eq. (1) see Paper I. 

The evolution of the chemical abundance distribution 
caused by diffusion processes during the whole white dwarf 
regime has been taken into account in this work. Our time- 
dependent element diffusion treatment, based on the for- 




Semiconvective mixing — see Straniero et al. (2003) 
not been considered in our models. 
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Fig. 8. Same as Fig. [T] but for the case of the OV model- 
sequence. 

mulation for multicomponent gases presented by Burgers 
(1969), considers gravitational settling, chemical and ther- 
mal diffusion'^ for the nuclear species ^H, "^He, ^He, ^^C, 
and ^^O. In this way, we avoid the use of the trace element 
approximation usually invoked in most ZZ Ceti studies. 

As for overshooting, we have included time-dependent 
overshoot mixing during all pre-white dwarf evolutionary 
stages. We have considered exponentially decaying diffu- 
sive overshooting above and below any formally convective 
region, including the convective core (main sequence and 
central helium burning phases), the external convective 
envelope and the short-lived helium-flash convection zone 
which develops during the thermal pulses. Specifically, we 
have followed the formalism of Herwig (2000). 

The effect of crystallization on the pulsational pattern 
of ZZ Ceti stars constitutes the central point of this pa- 
per. For a white dwarf model characterized by a given Tcff 
value, stellar mass and chemical composition, LPCODE 
provides a formally self-consistent value of the crystal- 
lized mass fraction (Mc/M*). In our calculations, crys- 
tallization is assumed to occur when F > 180 (Ogata & 
Ichimaru 1987; Stringfellow, De Witt & Slattery 1990), 
where F = Z'^e^ /rk-QT is the ion coupling constant. In or- 
der to evaluate self-consistently the perturbations caused 
by a crystal/fluid interface on the eigenfunctions we shall 



Radiative levitation has been neglected. 
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Fig. 9. Internal ^^O chemical profiles in terms of the frac- 
tional mass for the OV sequence corresponding to various 
percentages of crystallization. Thick-solid line corresponds 
to the profile just before the onset of crystallization, and 
thin lines correspond to the chemical profiles resulting 
from chemical rehomogenization at increasing percentages 
of crystallized mass fraction, from 10 to 40 % in steps of 
10 %. The curve at 40 % of crystallized mass fraction cor- 
responds to Tes ~ 10000 K. Thick dashed line represents 
the oxygen profile by the time the chemical discontinuity 
left by core overshooting has been completely wiped out 
by chemical rehomogenization (when « 37% of the core is 
in a crystalHzed state). 



adopt special boundary conditions (see §3.1). In addition, 
the chemical redistribution due to phase separation has 
been taken into account following the procedure described 
in Montgomery et al. (1999) and Salaris et al. (1997). 
To assess the enhancement of oxygen in the crystallizing 
core, we have employed the phase diagram of Segretain & 
Chabrier (1993). Details concerning the algorithm used to 
compute the mixing processes due to crystallization will 
be described in §3.5 below. 

An important aspect of the present study is to explore 
the implications of the occurrence of overshooting during 
the pre-white dwarf evolution for the pulsational proper- 
ties of crystallized ZZ Ceti stars. To this end, we consider 
the evolutionary sequences analyzed in Paper I, that is se- 
quence NOV based on the complete evolution of a 7.5-Mq 
initial mass star in which overshooting was not consid- 
ered, and sequence OV based on the evolution of a 6-Mq 
star with overshooting (a solar-like initial composition has 
been adopted). For both sequences, the mass of the result- 
ing carbon/oxygen core is quite similar (« 0.94 Mq). This 
has enabled us to compare the oscillation properties of 
crystallized ZZ Ceti stars characterized by the same stel- 



lar mass but being the result of the evolution of progenitor 
stars with different initial masses. It is worth mentioning 
that the ZZ Ceti models we employed in our analysis are 
the result of the complete evolution of massive white dwarf 
progenitors from the zero-age main sequence through the 
thermally pulsing and mass loss phases to the white dwarf 
regime. In particular, sequence OV experiences the third 
dredge-up and hot bottom burning during the thermally 
pulsing phase on the asymptotic giant branch (AGB). For 
details concerning the evolutionary properties of our se- 
quences we refer the reader to Paper I. In what follows, 
we restrict ourselves to commenting on the main results 
for the inner chemical composition. 

The upper panels of Figs, ^and El show the chemi- 
cal profiles at the ZZ Ceti stage for sequences NOV and 
OV. As demonstrated in Paper I, except for the inner part 
of the core, where the diffusion time scale becomes much 
longer than the evolutionary time scale, element diffusion 
is so efficient that the resulting abundance distribution at 
the ZZ Ceti stage does not depend on the occurrence of 
overshooting in the convective envelope during the ther- 
mally pulsing AGB phase. In fact, the resulting external 
chemical profile at the ZZ Ceti stage is quite similar for 
both sequences. Remarkably enough, our models are char- 
acterized by a chemical interface in which helium, carbon 
and oxygen in non- negligible abundances coexist, an in- 
terface which, when the white dwarf reaches the ZZ Ceti 
stage, has extended appreciably as a result of chemical dif- 
fusion. Needless to say, the trace element approximation 
would not be in this case an appropriate approach to treat 
element diffusion. 

However, overshoot episodes occurring during the core 
burning phases, particularly during the helium core burn- 
ing, leave recognizable features in the inner chemical 
profile of a massive ZZ Ceti star, as it is documented 
by Fig. 13 That is, the innermost region of such stars 
keeps a record of the extra mixing experienced during the 
pre-white dwarf evolution. In particular, note the sharp 
variation of the carbon/oxygen profile at logg w —0.7 
((7 = 1 — Mr / M^,) . As shown in Paper I, this is responsible 
for noticeable structure in the period spacing diagrams. 
The shape of the chemical profile towards the central re- 
gion of our OV model-sequence is typical for situations 
in which additional mixing beyond what is predicted by 
Schwarzschild criterion for convective stability is allowed. 
The occurrence of such mixing episodes, particularly core 
overshooting and/or semiconvection, is suggested by both 
theoretical and observational evidence. In particular, ex- 
tra mixing episodes beyond the external border of the fully 
convective core that take place towards the end of central 
helium burning have a large influence on the carbon and 
oxygen distribution in the core of white dwarfs. Recently, 
Straniero et al. (2003) have presented a detailed analysis 
of the inner chemical abundance in a 3-Mq star model 
resulting from different extra mixing processes during the 
late stage of core helium burning phase. In particular, they 
conclude that models which incorporate semiconvection 
or a moderate mechanical overshoot applied to core and 
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Fig. 10. Evolution of the periods for i = 2 modes as a 
function of the effective temperature, corresponding to 
the NOV model-sequence. The full effects of crystalliza- 
tion have been considered both in assessment the final 
shape of the internal chemical profile and in computing 
the eigenmodes. 

convective shells, predict a sharp variation of the chemi- 
cal composition in the carbon/oxygen core, in agreement 
with our results for OV sequence. 

3. Pulsational analysis 

In this section we shall explore in detail the consequences 
for the theoretical period spectrum of our white dwarf 
models due to the effects of crystallization. In all the cases, 
we shall analyze model-sequences NOV and OV in models 
with Toff decreasing from « 15000 to « 10000 K, which 
covers the temperature range for the observed ZZ Ceti in- 
stability strip. In showing our results we shall concentrate 
on £ = 2 f;-modes only; the results adopting other values 
of £ are qualitatively similar. 

Before discussing the pulsational results, we briefly de- 
scribe below the pulsational code and the treatment em- 
ployed to assess the various pulsational quantities. 

3.1. Pulsation code 

For computing adiabatic, nonradial y-modes of the ZZ 
Ceti models we employ the same pulsational code as in 
Corsico et al. (2001, 2002), with appropriate modifica- 
tions for handling the effects of crystallization on oscil- 
lation eigenmodes. Briefly, the code, which is coupled to 
the LPCODE evolutionary code, is based on the general 
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Fig. 11. Same of Fig. but for the OV model-sequence 



Newton-Raphson technique to solve the full set of lin- 
earized equations governing spheroidal"*, adiabatic, non- 
radial pulsations of spherically symmetric stars. The code 
provides the eigenperiod Pk (being k the radial overtone 
of mode) and the dimensionless eigenfunctions yi, • • ■ , 2/4 
— see Unno et al. (1989) for their definition. Useful pulsa- 
tion quantities, such as the oscillation kinetic energy, the 
weight function, and the variational period — see Kawaler 
et al. (1985) — are also provided by our pulsational code 
for each computed eigenmode. Finally, the asymptotic pe- 
riod spacing is computed as in Tassoul et al. (1990). 

The boundary conditions arc those given by Osaki & 
Hansen (1973) and the normalization condition adopted 
at the stellar surface is yi = 1. As for the inner boundary 
conditions of our crystallizing models, we have adopted 
the "hard-sphere" boundary condition, which has been 
shown by Montgomery & Winget (1999) to be a realistic 
representation. In fact, the amplitude of eigenfunctions of 
(7-modes is drastically reduced below the solid/liquid in- 
terface due to the non-zero shear modulus of the solid, as 
compared with the amplitude in the fluid region — see 
figure 4 of Montgomery & Winget (1999). 



^ We have not considered torsional modes, since these modes 
are characterized by very short periods — up to 20 s; see 
Montgomery & Winget (1999) — and are not observed in ZZ 
Ceti stars. 
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Fig. 12. The forward period spacing in terms of the pe- 
riods for £ — 2 modes corresponding to the NOV model- 
sequence analyzed in Fig. IIUI Each panel shows the pre- 
dictions for different percentages of the crystallized mass 
fraction. Here chemical rehomogenization induced by crys- 
tallization has been considered. 



Specifically, the hard-sphere boundary condition at the 
radial shell corresponding to the outward-moving crystal- 
lization front (rc = r{Mc)) reads 

yi = 0, 

j/2 = arbitrary, (2) 
^ ys - 2/4 = 0. 

Here, yi and y2 represent the radial and tangential dis- 
placements, respectively, and and y^ are the Eulerian 
perturbation of the gravitational potential and its deriva- 
tive. Note that the last condition is the same as for the 
normal case in which the core is in a fluid state and the 
boundary condition is applied at the stellar centre — see 
Appendix B of Montgomery & Winget (1999). 

3.2. The Brunt-Vaisala frequency 

The procedure we follow to assess the run of the Brunt- 
Vaisala frequency (N) is that proposed by Brassard et al. 
(1991b). This numerical treatment accounts explicitly for 
the contribution to N from any change in composition in 
the white dwarf models by means of the Ledoux term B. 
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This is an important aspect in connection with the phe- 
nomenon of mode trapping and confinement (Brassard et 
al. 1992a; Corsico et al. 2002). In the lower panels of Figs, 
n and 12 we show the run of N'^ in terms of the outer 
mass fraction for white dwarf models at Toff ~ 11800 K 
corresponding to sequences NOV and OV, respectively. 
Element diffusion strongly smooths out the external chem- 
ical profile emerging from the thermally pulsing AGB 
phase to such an extent that the Brunt-Vaisala frequency 
in the outer layers exhibits very smooth local features. 
Note that because of the markedly distinct shape of the 
chemical distribution in the innermost region of the model 
with core overshooting, the profile of N for such model ex- 
hibits a clearly peaked feature at logg ^ —0.7 (Fig. 
This feature is responsible for the presence of pronounced 
minima in the period spacing distribution (see figure 15 
of Paper I). The model corresponding to sequence NOV 
also shows an innermost peak in N (Fig. However, 
because it is located at a deeper layer and has a smaller 
amplitude than in the OV sequence, it produces a much 
smaller amount of mode trapping. 

3.3. Results without crystallization 

In this section we shall describe pulsational results in 
which we have neglected crystallization in the computa- 
tion of the pulsational spectrum. In Fig. 13 we show the 
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evolution of the i — 2 periods in terms of Toff correspond- 
ing to the NOV sequence of models. Clearly, the periods 
increase monotonically with decreasing Teg, particularly 
in the case of modes with high radial order k. This ef- 
fect is explained on the basis that the Brunt- Vaisala fre- 
quency at the core decreases as a result of increasing de- 
generacy in the core. Note that the mode with fc = 1 is 
rather insensitive to the white dwarf cooling, though its 
period slightly grows with decreasing Teff. We note that 
the curves are very smooth in all the T^s range consid- 
ered, and only in a few cases is there a slight approach 
between periods of adjacent modes. In Fig. ^we depict 
the evolution of periods in the case of models correspond- 
ing to OV sequence. In this case the overall trend of the 
periods is similar to that found in NOV models, but they 
show instead clear signals of mode bumping and avoided 
crossing^. For instance, a series of mode bumping is ob- 
served between modes with fc = 17 and fc = 18, fc = 18 
and A: = 19, A: = 19 and fc = 20, and fc = 20 and fc = 21 
(thick lines in Fig. We also find that, apart from the 
periods, the kinetic energy and period spacing minima are 
also exchanged between these pairs of modes after mode 
bumping has taken place. A close inspection of the eigen- 
functions and weight functions demonstrates that modes 
experiencing avoided crossings are core trapped ones, that 
is, modes with relatively large amplitudes in the region 
bounded by the centre of the model and the location of 
the step in the oxygen profile left by overshooting (see 
Paper I). Similar mode bumping/ avoided crossing features 
in white dwarf evolutionary models have been reported by 
Wood & Winget (1988) and Brassard et al. (1991a), al- 
though in such calculations modes experiencing avoided 
crossings were identified as modes trapped in the outer 
hydrogen envelope. In the context of massive white dwarf 
models, Montgomery & Winget (1999) have found avoided 
crossing phenomena when the crystallized mass fraction 
of their models is varied in a continuous fashion and the 
other model parameters (such as Tos) are held fixed. 

3.4. Results considering crystallization without 
chemical rehomogenization 

In this section, we describe our results when account is 
made of the presence of a solid core in the computation 
of the pulsational eigenspectra. However, we neglect here 
the effect of chemical rehomogenization due to phase sep- 
aration, that is, the core chemical profile remains fixed 
as the white dwarf crystallizes. We stress that we take 
the crystallized mass fraction not as a free parameter but 
as given by evolutionary computations. Figs. and 
depict the resulting evolution of periods in terms of Tcff 
for sequences NOV and OV, respectively. Note that in 
both cases the onset of crystallization is reflected in the 
pulsational spectrum as a rather abrupt increase in peri- 

^ When a pair of modes experiences avoided crossing, the 
modes exchange their intrinsic properties after mode bumping; 
see Aizenman et al. (1977). 



ods for modes with high- and intermediate radial orders, 
whereas lower overtones are rather insensitive to the pres- 
ence of the solid core. This fact can be explained on the 
basis that for (^-modes, when fc ^ 1 (asymptotic limit) 
the periods are given by Pfc oc fc (/J'" Ndr/r)-'^ (Tassoul 
et al. 1990), where ri and r2 define the propagation re- 
gion of modes. When crystallization is considered, the in- 
ternal boundary ri(Mc) moves outward, so the integral 
f^^ Ndr/r decreases, and consequently the periods (and 
the period spacings) increase. An important feature dis- 
played by Figs. |5l and is that, when crystallization 
begins, the mode bumping/ avoided crossing phenomena 
propagate to longer periods. Indeed, the outward-moving 
crystallization front reinforces the process of mode bump- 
ing/avoided crossing, which are also present before the 
onset of crystallization. This effect is particularly notice- 
able in the case of OV models. In Figs. Qand |H|we show 
the period spacing diagrams corresponding to various per- 
centages of the crystallized mass fraction for NOV and OV 
model sequences, respectively. Note that in the case of OV 
models the strong minima in the period spacing are grad- 
ually displaced to longer periods as crystallization pro- 
ceeds. At the same time they move away from each other. 
Note that, at T^s ~ 10000 K (40% of crystallized mass 
fraction) only one minimum at about 550 s remains. The 
period spacing diagrams exhibited by models with core 
overshooting are completely different from those of models 
in which overshooting is neglected, as can be inferred from 
Fig. In fact, although the period spacing diagrams are 
strongly modified by crystallization, their overall struc- 
ture looks very different between both sequences. In par- 
ticular, sequence NOV lacks the strong minima character- 
izing sequence OV, features which, as concluded in Paper 
I, could be eventually used to place constraints on the oc- 
currence of core overshooting. So, the conclusions arrived 
at in Paper I concerning the role of core overshooting in 
massive ZZ Ceti stars remain valid even when the effect 
of a solid core on the pulsations is included. However, this 
comes with the caveat that if the star is more massive than 
this model, then it could have a large enough crystallized 
mass fraction (^ 80 %) so that the chemical step left by 
overshooting would be engulfed by the crystallizing core, 
erasing this feature. In this case we should expect a pulsa- 
tional pattern with no appreciable minima in the period 
spacing. As we will see in the following section, this situa- 
tion can be reached for smaller amounts of crystallization 
if the effects of phase separation and mixing are taken into 
account. 

3.5. Results considering crystallization and chemical 
rehomogenization 

As recently shown, the occurrence of core overshooting 
strongly modifies the appearance of the period diagrams 
even in the presence of a solid core. Now, theoretical 
evidence strongly suggests that, if the white dwarf core 
is composed initially of a mixture of carbon and oxy- 
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Fig. 14. Period spacing vs. period for effective temperature values between w 15000 and 10000 K corresponding to 
the NOV model-sequence. In the interests of clarity, each curve has been shifted upward starting from the lowest one, 
in increments of 1 sec. For details, see the text. 



gen, then the crystallized region will have an enhanced 
abundance of oxygen compared to that in the original 
fluid state (Ichimaru et al. 1988; Segretain & Chabrier 
1993, Salaris et al. 1997). On the other hand, the fluid 
regions overlying the crystallized oxygen-enhanced layers 
will have a higher content of carbon. This region will be- 
come Ray leigh- Taylor unstable, since carbon is slightly 
less dense than oxygen. This instability leads to a reho- 
mogenization of the chemical profile, and, as a result, the 
chemical abundances after crystallization can differ sub- 
stantially when compared with the initial ones. This is 
an important point because any feature present in the 
chemical profile could be potentially wiped out by this 
rehomogenization process, even in fluid layers located far 
away from the crystallization front. 

To compute the above-mentioned chemical rehomog- 
enization, we have employed the same algorithm as in 
Montgomery et al. (1999) and Salaris et al. (1997). To 
derive the oxygen enhancement when a given layer crys- 
tallizes, we adopt the Segretain & Chabrier (1993) spindle- 
type phase diagram for a carbon/oxygen mixture. In our 
scheme for mixing, we first consider a crystallized, oxygen- 
enhanced layer and then we ask if the innermost fluid shell 
has a higher carbon content than the overlying one. If it 
does, then we mix both fluid layers and perform the same 



comparison with the next layer further out. In this way 
we scan outward through the fluid and the process stops 
when further mixing no longer decreases the carbon con- 
tent of the fluid between this point and the crystallization 
boundary. When the crystallization front moves outward 
due to cooling, the aforementioned procedure is repeated. 
In Fig. 1^ we show the oxygen chemical profile in terms 
of the fractional mass for the OV model sequence corre- 
sponding to various degrees of crystallization. The initial, 
uncrystallizcd profile is depicted with a thick-solid line, 
and the resulting profiles after the rehomogenization has 
taken place are shown with thin lines for 10, 20, 30 and 
40 % crystallized mass fractions. Note that the shape of 
the chemical profile is strongly modified, even in regions 
far from the crystallization front. As a clear example of 
this, we show with thick-dashed line the oxygen profile at 
the moment that the pronounced feature at Mr « 0.80Af* 
(left by core overshooting) is completely eliminated. This 
occurs when the location of the crystallization boundary 
reaches « 0.37Af*. 

We wish now to investigate how such chemical reho- 
mogenization affects the pulsational spectrum of the white 
dwarf models. Figs. 1101 and 1111 document the evolution 
of ^ = 2 periods in terms of Tcs corresponding to cases 
NOV and OV sequences, respectively. In the case of the 
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Fig. 15. Same as Fig. EI but for the OV model-sequence. 



NOV sequence, the rehomogenization process induced by 
crystallization clearly modifies the internal chemical pro- 
file (it smooths out the profile at logq w —0.3) to such a 
degree that any signal of mode bumping/avoided cross- 
ing is almost absent, in contrast with the situation in 
which such mixing has been neglected (see Fig. ISJ- In 
the case of OV models, the theoretical spectrum of peri- 
ods barely changes in response to the rehomogenization 
(compare Figs. Eland That is, clear features of mode 
bumping/avoided crossing survive even during the stages 
in which the pronounced step in the oxygen profile is con- 
tinuously diminished by chemical rehomogenization (see 
Fig. However, a closer inspection of Fig. 1111 reveals 
that from Mc/Af* « 0.37 (Toff « 10400 K) these features 
virtually disappear, and the period spacing becomes al- 
most uniform. This effect can be appreciated more clearly 
from Fig. 1131 in which we have plotted period spacing 
diagrams for several degrees of crystallization. Note that 
the pronounced minima characterizing the period spacing 
distribution at periods ^ 760 and ^ 1120 s (for 30 % crys- 
tallization) are absent when the crystallized mass fraction 
is 35%. In fact, when the chemical discontinuity has been 
completely wiped out by chemical rehomogenization, the 
predicted period spacing distribution becomes very flat, 
particularly for periods longer than ^ 300 sec. 

Figs. 1 121 and El represent some of the main results of 
the present work. A comparison of these figures clearly in- 



dicates that the period spacing distribution for sequences 
NOV and OV is markedly different. In particular, we note 
that in the case of the OV model with Ale = 0.35M*, after 
the overshoot-induced step has been almost removed from 
the chemical profile, the period spacing distribution is no- 
tably more uniform as compared with the NOV model 
characterized by the same percentage of crystallization. 
This is because in regions farther out from the centre 
(logg ^ —1) the OV model is actually smoother than 
the NOV model, so its period spectrum shows less mode 
trapping. 

Finally, in Figs. 1141 and 1151 we show the forward pe- 
riod spacing vs period for 15000 K ^ Toff <; 10000 K cor- 
responding to the NOV and OV model-sequences, respec- 
tively. For the sake of clarity, each curve has been shifted 
upward starting from the lowest one, in increments of 1 
sec. Labels of the solid curves indicate the percentage of 
crystallization and the corresponding Toff of the models. 
The dashed curve corresponds to the moment at which the 
relevant structure in the oxygen profile (at logq « —0.3 in 
the case of the NOV sequence, and at \ogq « —0.7 in the 
case of the OV sequence) disappears due to chemical reho- 
mogenization. From these plots one can see the evolution 
of the trapped modes quite clearly, particularly when crys- 
tallization begins. Note from Fig. El the decrease in am- 
plitude of mode trapping as the degree of crystallization 
increases. This figure also emphasizes the gradual change 
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Fig. 16. Top panel: the bumps/non-smooth features in 
the Brunt- Vaisala frequency found in our numerical mod- 
els (dashed curve) and the assumed bumps which we use 
in our simple analytical model (solid curve). In the lower 
panels, we illustrate how these features can act to trap 
modes. For instance, the narrowness of the $ = 0.16 fea- 
ture, compared to the spatial wavelength of these modes, 
causes it to induce amplitude changes in the two eigen- 
modes pictured. 



experienced by the period spacing when the step left by 
overshooting is removed by chemical rehomogenization. 



4. Analytical model of mode trapping 

In this section we examine the mode trapping in our mod- 
els by comparing it to the results of a simple analytical 
model. This analytical model, which will be described in 
more detail elsewhere, describes the high-overtone limit of 
the pulsation equations as derived by Deubner & Gough 
(1984) and Gough (1993). While formally equivalent to the 
standard linear adiabatic wave equation in the Cowling 
approximation — see, e.g., Cox (1980) — this form of the 
equations reduces the problem to a single second-order 
differential equation resembling the oscillations of a non- 
uniform string. In fact, using this equation as a convenient 
starting point, Montgomery, Metcalfe, & Winget (2003) 
showed that a symmetry inherent in the string problem is 
still present in stellar pulsations. This lead to the discovery 
of the core/envelope symmetry and the subsequent real- 
ization that the chemical composition gradients in the core 
may be even more effective in producing mode trapping 
than those in the envelope. Here we use an extended ver- 
sion of this model in which the perturbations to the back- 
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Fig. 17. Comparison of the mode trapping of the analyt- 
ical model (filled circles, solid lines) to that of the full 
numerical results shown in Fig. 13 (open circles, dotted 
lines). Each panel is labelled with the crystallized mass 
fraction to which it corresponds. 

ground state need not be small, although our solutions are 
still formally valid only in the high-overtone limit. 

In Fig. El we illustrate how features/bumps in the 
Brunt- Vaisala frequency result in mode trapping in the 
eigenf unctions. As in Montgomery et al. (2003), we use 
the "normahzed buoyancy radius", defined by 

<i>(r) ^ - 



as the radial variable, since with this choice the wavelength 
of the eigenfunctions appears to be constant, making the 
correspondence with the vibrating string more apparent 
(<& equals and 1 correspond to the centre and surface 
of the model, respectively). In addition, the eigenfunction 
vl* « {r'^ /gp^^^) 5p is defined in such a way that it has con- 
stant amplitude except in regions where the background 
state varies rapidly, such as the bumps shown in the top 
panel of Fig. El As can be seen, the "bump" at $ ^ 0.16 
is narrower than the wavelengths of either the A; = 5 or 
A; = 30 mode, so both modes perceive it essentially as a 
delta function and therefore experience significant mode 
trapping (amplitude change) in propagating through it. In 
contrast, the other two "bumps" at $ 0.3 and $ ~ 0.6 
are much wider. Given that the wavelength of the fc = 30 
mode is smaller than the width of these bumps, it expe- 
riences essentially no mode trapping. On the other hand, 
the wavelength of the k = b mode is somewhat larger than 
the width of these bumps, so it still experiences a modest 
amount of mode trapping (amplitude change) due to these 
features. 
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In the top panel of Fig. ^| we have chosen the height, 
width, and position of the bumps in the analytical model 
(the solid curve) to mimic those found in a full evolution- 
ary model (the dashed curve). In Fig.^lwe show the mode 
trapping found in our analytical model (filled circles, solid 
lines) and that found in the OV model shown in Fig. 13 
(open circles, dotted lines), for values of the crystallized 
mass fraction between % and 40 %. Clearly, the over- 
all correspondence is quite good, except for low periods, 
where our asymptotic model is anyway not applicable. In 
particular, the asymptotic model is able to reproduce i) 
the shape of the trapping cycle (sharp minima and flat 
maxima), ii) the decrease in amplitude of mode trapping 
for high periods/overtones, iii) the increase in the period 
of the trapping cycle and the movement of mode trapping 
features to higher periods as the degree of crystallization 
is increased, and iv) the absence of trapping features for 
the 40 % crystallized case. 

The explanation of these effects is straightforward. 
First, the overall good match is due to the fact that the 
modes shown actually have high overtone numbers and 
that the mode trapping is dominated by the single feature 
in SN/N at $ 0.16. Second, the decrease in trapping 
with increasing period occurs because these modes have 
shorter spatial wavelengths, and when these wavelengths 
become comparable to the width of the features in 6N/N 
the modes cease to feel the effects of a sudden transition. 
Third, the reason that the number of modes per mode 
trapping cycle increases is that crystallization effectively 
moves the $ = point to the right, making the bump at 
$ ~ 0.16 closer to the inner turning point. The shortest 
possible trapping cycle of 2 modes is achieved for a feature 
in the middle at $ ~ 0.5, so moving a bump closer to an 
edge has the effect of making the trapping cycle longer. 
Finally, the lack of mode trapping in the 40 % crystallized 
case is due to the removal of the $ ~ 0.16 bump due to 
crystallization, phase separation, and the subsequent re- 
homogenization: the remaining bumps cause virtually no 
mode trapping. 

As a final note we point out that the core/envelope 
symmetry found by Montgomery et al. (2003) is not an 
issue for these models, at least for modes with periods 
greater than 300 sec. This is because, as is conclusively 
demonstrated in Fig. 17, the outer chemical transition 
zones (C/He and He/H) are found to be so smooth that 
they produce almost no mode trapping; thus, any mode 
trapping must be a result of structure in the deeper re- 
gions. 

5. Discussion and conclusions 

Convective overshooting is a longstanding problem in the 
theory of stellar structure and evolution. It is well known 
on theoretical grounds that during many stages in their 
lives stars experience overshoot episodes, that is, partial 
mixing beyond the formally convective boundaries as pre- 
dicted by the Schwarzschild criterion of convective sta- 
bility (Zahn 1991; Canuto 1992; Freytag et al. 1996; see 



also Renzini 1987). In particular, core overshooting taking 
place during central burning is an important issue because 
it has significant effects on the stellar structure and evo- 
lution. Over the years, considerable observational effort 
has been devoted to demonstrating the occurrence of core 
overshooting. Indeed, confrontation of stellar models with 
a wide variety of observational data suggests that convec- 
tive overshoot takes place in real stars — see Stothers & 
Chin (1992), Alongi et al. (1993), Kozhurina-Platais et al. 
(1997), Herwig et al. (1997), von Hippel & Gilmore (2000), 
among others. 

However, most of the evidence about the occurrence 
of core overshooting relies primarily on observational data 
from the very outer layers of stars from where radiation 
emerges. A more promising and direct way of placing con- 
straints on the physical processes occurring in the very 
deep interior of stars is by means of the study of their 
pulsational properties. Pulsating white dwarfs are partic- 
ularly important in this regard. In fact, white dwarfs con- 
stitute the end product of stellar evolution for the vast 
majority of stars, and the study of their oscillation spec- 
trum through asteroseismological techniques has become 
a powerful tool for probing the otherwise inaccessible in- 
ner regions of these stars (Bradley 1998b; Metcalfe et al. 
2002; Metcalfe 2003). White dwarf asteroseismology has 
also opened the door to peer into the physical processes 
that lead to the formation of these stars. 

In this work — and also in Paper I — we have argued 
that the occurrence of core overshoot episodes during core 
helium burning leaves strong imprints on the theoretical 
period spectrum of massive ZZ Ceti stars, features which 
could in principle be used for providing strong constraints 
on the occurrence of such episodes. On the other hand, 
we have demonstrated that the chemical rehomogeniza- 
tion induced by phase separation gives rise to a feature- 
less period spectrum when the crystallized mass fraction 
is larger than « 30%. One of the aims of the present paper 
was to place constraints on the stellar mass and T^ff values 
at which we should expect a pulsational pattern without 
any signature of core overshoot. We have seen that in our 
sequence OV the chemical discontinuity left by core over- 
shooting is wiped out by chemical rehomogenization for 
a solid core larger than 37%, which for a 0. 94- A/q white 
dwarf takes place at Tcs — 10400 K. The effective tem- 
perature at which this occurs is strongly dependent on 
the stellar mass. For instance, for 1.00-, 1.03-, and 1.05- 
Mq white dwarf models, the traces of core overshooting are 
removed by chemical rehomogenization at Tcff — 11700, 
12350, and 12800 K, respectively 

The white dwarf star BPM 37093, the most massive 
ZZ Ceti known to date, is particularly noteworthy in this 
regard. Using lUE spectra, Koester & AUard (2000) derive 
Teff= 11520 K, log 5 = 8.67 and = 1.03 Mq. Other 
recent determinations place the mass of BPM 37093 at 
1.00 M0 (Bergeron et al. 2001) and 1.10 M© (Bergeron 
et al. 2004; Fontaine et al. 2003). Even given these un- 
certainties, it seems likely that its mass lies in the range 
1.0-1.1 Mq. With regard to its pulsation properties, BPM 
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37093 exhibits g-mode pulsations with periods in the 
range 500 — 660 s (Kanaan et al. 2000). The periods are 
very accurately determined, but as far as wc arc aware, the 
£ values are not securely identified^. On general grounds, 
it is possible to show that the observed modes cannot all 
be £ = 1, so the simplest assumption has been that they 
are ^ = 2 modes'^ (Nitta et al. 2000, Kanaan et al. 2000, 
Montgomery & Winget 1999). Under this assumption, it 
is possible to construct an "observed" period spacing di- 
agram, which is characterized by an average period spac- 
ing of 17.125 s with deviations due to mode trapping of 
» ±8 s. 

Given the uncertainties in the mass and Tqs of BPM 
37093 and also in its mode identification, we cannot make 
a definite statement concerning cither whether core over- 
shooting took place in BPM 37093 or what degree of 
crystallization it possesses. However, we can place con- 
straints on the different possibilities. In fact, we can 
characterize BPM 37093 by a stellar mass in the range 
1.00 Mq ^ M* ^ 1.10 Mq and an effective temperature 
in the interval 11000 < Tcff <, 12000 K, given the present 
uncertainties in the derivation of these quantities. In the 
case of 1.00 Mq our treatment of crystallization predicts 
that the model should be 37% crystallized at « 11700 K. 
Thus, its observed frequency spectrum could reflect mode 
trapping effects due to prior evolutionary episodes of core 
overshoot. On the other hand, if BPM 37093 has a mass 
of PS 1.10 Mq, then there should be virtually no mode 
trapping features in its period spectrum. In this case, we 
cannot say anything about the occurrence of core over- 
shoot during pre-white dwarf evolution. A possible expla- 
nation of the observed period spectrum would be that the 
observed modes are not all £ = 2, m=0 modes, but instead 
a mixture oi £ = 1 and £ = 2 modes (Metcalfe et al. 2004, 
Nitta et al. 2000). 

6. Concluding remarks 

This paper is aimed at specifically assessing the feasibil- 
ity of employing white dwarf asteroseismology to demon- 
strate the occurrence of core overshooting during core 
helium burning. The investigation is focused on mas- 
sive, intermediate-mass stars, the complete evolution of 
which has been followed from the zero-age main sequence 
through the thermally pulsing phase on the asymptotic gi- 
ant branch (AGB) to the white dwarf regime. This allows 
us to obtain ZZ Ceti models with a physically sound in- 
ternal structure consistent with the predictions of stellar 
evolution theory. Our work shows that for studying the 
effects of core overshooting, pulsating white dwarfs with 

^ By applying the £- value identification method of Robinson 
et al. (1995), Nitta et al. (2000) have discarded the possibility 
that the main periods observed in BPM 37093 are £ = 3 modes, 
but they cannot give definite identifications regarding ^ = 1 or 
£ = 2 values. 

^ If all the modes are £ = 1 then the average period spacing 
would imply a mass much larger than the values found from 
the spectroscopic determinations. 



masses of ~ 0.90 Mq would be ideal since the potential for 
crystallization will not be a factor in erasing the signature 
of these overshooting events. Conversely, for studying the 
effects of crystallization, masses ^ 1.0 M© are best since 
most of the core structure will be erased by crystalliza- 
tion and the chemical mixing associated with phase sepa- 
ration, making it easier to isolate and study the effect of 
crystallization on its own. In future studies we will focus 
on models in these mass ranges. 
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